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Abstract — In this paper, two low-complexity adaptive step size 
algorithms are investigated for blind adaptive beamforming. Both of 
them are used in a stochastic gradient (SG) algorithm, which employs 
the constrained constant modulus (CCM) criterion as the design 
approach. A brief analysis is given for illustrating their properties. 
Simulations are performed to compare the performances of the novel 
algorithms with other well-known methods. Results indicate that the 
proposed algorithms achieve superior performance, better convergence 
behavior and lower computational complexity in both stationary and 
non-stationary environments. 

Index Terms-Blind adaptive beamforming techniques, constrained 
constant modulus (CCM), modified adaptive step size (MASS), time 
averaging adaptive step size(TAASS). 

I. Introduction 

In recent years, many adaptive filtering algorithms have been 
used for beamforming, which is a promising and widely investigated 
technology for rejecting interference and improving the performance 
of high capacity mobile communications systems |1|. Many methods 
have been presented in different communication systems [2|-j4|- In 
contrast to fixed beamforming techniques, an adaptive beamformer 
has the ability of rejecting interference and pointing its mainbeam 
in the desired direction with the change of scenarios. Blind adaptive 
beamforming, which is intended to form the array direction response 
without knowing users' information beforehand, is an important topic 
that deals with interference cancellation, tracking improvement and 
complexity reduction. 

The blind adaptive SG method, which is commonly employed in 
the blind adaptive beamforming area, is a well-known technique for 
solving optimization problems with different criteria, e.g., minimum 
mean squared error (MMSE) j3], minimum variance (MV) [6] and 
constant modulus (CM) |7|. The MV algorithm is a computational 
efficient approach to estimate the input covariance matrix. The results 
in [6 1 prove that the MV criterion leads to a solution identical to that 
obtained from the minimization of the mean squared error (MSE). 
The CM algorithm for beamforming exploits the low modulus fluc- 
tuation exhibited by communications signals using constant modulus 
constellations to extract them from the array input. It is well known 
that the performance of the CM method is superior to that of the 
MV. A disadvantage of both two methods is that they are sensitive 
to the step size. The small value of the step size will lead to slow 
convergence rate, whereas a large one will lead to high misadjustment 
or even instability. Besides, the CM cost function may have local 
minima, and CM receivers do not have closed-form solutions. Xu and 
Liu (8 1 developed a SG algorithm on the basis of the CCM technique 
to sort out the local minimum problem and obtain the global minima. 
But the problem incurred by the step size cannot be solved properly. 

For accelerating the convergence, recursive least squares (RLS) 
algorithms were introduced by Xu and Tsatsanis using the constrained 
minimum variance (CMV) criterion [9] and then developed by de 
Lamare and Sampaio-Neto with the CCM approach [10]. The latter, 
which improves the performance significantly, optimizes a quadratic 
cost function based on the CM criterion subject to linear constraints 
for the array weight adaptation. Combining with the constrained 
condition, this method reaches an optimal solution. Nevertheless, 



the RLS based beamformers cannot avoid complicated computations 
caused by the required correlation matrix inversion. 

Comparing SG algorithms, which represent simple and low- 
complexity solutions but subject to slow convergence, with RLS 
methods, which possess fast convergence but high computational 
load, it is preferable to adopt SG beamformers due to complexity and 
cost issues. For this reason, the improvement of blind SG techniques 
is an important topic and has been investigated for several decades. 
In this research area, the works in |(6] and J9] employ standard SG 
methods with fixed step size (FSS) that are not efficient with respect 
to convergence and steady-state performance. The performance of 
the beamformer is strongly dependent on the choice of the step 
size [12]. It reflects a tradeoff between misadjustment and the 
convergence rate. Actually, the communication systems are non- 
stationary environments, which make it very difficult to predetermine 
the step size. It is necessary to make the step size track the change of 
the system automatically and so obtain good convergence behavior. 
Previous researches have focused on this aim and some good results 
have been reported. The adaptive step size (ASS) mechanism was 
employed in both the MV (T3), [TJ] and the CM (T5) criteria for 
improving the performance. Because of requiring an additional update 
equation for the gradient of the weight vector with respect to the step 
size, which increases the extra computational load, the applications of 
these algorithms are limited. The authors of [16] propose two novel 
variable step size mechanisms for MV algorithms. The simulation 
experiments show that the new mechanisms are superior to previously 
reported methods and have a reduced complexity. 

This paper proposes two blind CCM beamformers based on two 
novel adaptive step size mechanisms. The origins of these mecha- 
nisms can be traced back to the works of [17] and [18| where low- 
complexity adaptive step size mechanisms were developed for LMS 
algorithms. In contrast to [ 161- 11181 . the mechanisms here are designed 
for CCM algorithms, since it is well-known that they outperform the 
CMV algorithms [8|. Because of this reason, the simulations here just 
compare the algorithms related to the CCM criterion. The additional 
number of operations of the proposed algorithms does not depend on 
the number of sensor elements. In addition, the results are presented 
for stationary and non-stationary environments, proving that the new 
mechanisms could reach better performance and faster convergence 
behavior than those of previous methods. 

The remaining of this paper is organized as follows. In the next 
section, we present a system model for smart antennas. Based on 
this model, the blind adaptive CCM beamformer design using the 
SG method is presented in Section III. In Section IV, the proposed 
adaptive step size mechanisms are derived. Simulation results are 
provided in Section V, and conclusions are drawn in Section VI. 

II SYSTEM MODEL 

In order to describe the system structure, let us make two simpli- 
fying assumptions for the transmitter and receiver models [ 19]. First, 
the propagating signals are assumed to be produced by point sources; 
that is, the size of the source is small with respect to the distance 
between the source and the sensors that measure the signal. Second, 
the sources are assumed to be in the "far field," namely, at a large 



distance from the sensor array, so that the spherically propagating 
wave can be reasonably approximated with a plane wave. Besides, 
we assume a lossless, nondispersive propagation medium, i.e., a 
medium that does not attenuate the propagating signal further and 
the propagation speed is uniform so that the waves travel smoothly. 
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Fig. 1. Adaptive beamforming structure for ULA. 

Let us consider the adaptive beamforming scheme in Fig. [TJ and 
suppose that q narrowband signals impinge on the uniform linear 
array (ULA) of m (q < m) sensor elements from the sources with 
unknown directions of arrival (DOAs) 9o,. ■ ■ ,9 q -i- The ith snapshot's 
vector of sensor array outputs can be modeled as [20] 



x(i) 



A(0)s(i)+n(i), i = l,...,N (1) 

qxl is the vector of the unknown 

,a(6 q -i)] G C mxq is the 



where = [9 , . . . , 9 q -i] 1 € C 
signal DOAs, A(9) = [a(9o), ■ ■ ■ , a{o q -i 
complex matrix composed of the signal direction vectors a(9k) = 

[1 e ~ 27rj Tt C ° aek e -2Tj(™-l)^cosS fc jT £ Q-mxl Q. _ 

0, ...,q — 1), where A c is the wavelength and d = A c /2 is the 
inter-element distance of the ULA, s(i) G lZ qxl is the real value 
vector of the source data, n(i) G C mx is the complex vector of 
white sensor noise, which is assumed to be a zero-mean spatially 
and white Gaussian process, N is the number of snapshots, and (-) T 
stands for the transpose. The output of a narrowband beamformer is 
given by 



y(i) = w(i) x(i) 

IT £ Qmxl 



(2) 



where w(i) = hvi(i), . . . ,ui m (i)] J G Q-rnxi j s mg conl pj ex we ight 
vector, and (-) H stands for the Hermitian transpose. 

III. Blind Adaptive CCM Algorithms 

The purpose of SG algorithms is to get an acceptable output per- 
formance and reduce the complexity load by avoiding the correlation 
matrix estimation and inversion. We describe the CCM algorithm on 
the basis of the SG method, which is called CCM-SG. 

For the CCM-SG algorithm, we consider the cost function as the 
expected deviation of the squared modulus of the array output to a 
constant, say 1. The CCM cost function is simply a positive measure 
of the average amount that the beamformer output deviates from the 
unit modulus condition [4|. By using the constraint condition, the 
cost function of CCM-SG can be expressed as 



Jck 



|y(*)| 2 -l) 2 , i = l, 



subject to w (i)a(6*o) = 1 



(3) 



where 0(6*0) denotes the steering vector of the desired signal. The 
constrained optimization means that the technique minimizes the 
contribution of the undesired interferences while maintaining the gain 
along the look direction to be constant. 

The SG blind adaptive algorithm optimizes the Lagrangian cost 
function described by 



Lccm = (\y(t)f - if + \(w (i)a(9 ) - 1) 



(4) 



where A is a scalar Lagrange multiplier. The solution can be obtained 
by setting the gradient terms of (|4j with respect to io(i) equal to zero 
and using the constraint. Thus, we obtain 



w(i + 1) = w(i) - n{\y{i)\ 2 - l)y*(i) 
■ [x(i) - a H (9 )x(i)a(9 )] 



(5) 



where /1 here is the step size, which is a fixed value for FSS and a 
variable value for ASS and * denotes complex conjugate. 

Because of the shortcomings introduced before for both FSS 
and ASS algorithms, it is necessary to develop other methods for 
improving the performance of SG method. 

IV. Proposed Adaptive Step Size Mechanisms 

In this section, two novel adaptive step size methods are described 
for adjusting the step size following the change of the communi- 
cation system. The step size adjustment is controlled by the square 
prediction error, which means that a large error will cause the step 
size to increase for providing fast tracking while a small error will 
result in a decrease of step size to yield smaller misadjustment. The 
computational complexity is not a problem in these mechanisms. 

A. Modified Adaptive Step Size (MASS) Mechanism 

The first proposed algorithm based on the MASS mechanism 
employs the prediction error and uses the update rule 



^ + l) = a^i) + 1 (\y(i)\ 2 -if 



(6) 



where < a < 1, 7>0 and y(i) is the same as that in (|2j- The 
rationale for the MASS is that for large prediction error the algorithm 
will make the step size increase to track the change of the system 
whereas a small error will result in a decrease of the step size. The 
parameter 7 is an independent variable for controlling the prediction 
error and scaling it at different levels. It is worth pointing out that 
the step size fj,(i + 1) should be limited in a range as follows 



H(i + 1) 




if fi(i + 1) > A* n 
if n(i + 1) < n„ 
otherwise 



(7) 



where < fimin < Umax- The constant fi m in is chosen as a com- 
promise between the satisfying level of steady-state misadjustment 
and the required minimum level of tracking ability while yL max is 
normally selected close to the point of instability of the algorithm for 
providing the maximum convergence speed. The MASS is the result 
of several attempts to devise a simple and effective mechanism. 



TABLE I 

Additional Computational Complexity of Proposed 
Algorithms 



Proposed 
Algorithms 


Number of operations per snapshot 


Additions 


Multiplications 


MASS 


1 


3 


TAASS 


2 


6 



B. Time Averaging Adaptive Step Size (TAASS) Mechanism 

The second mechanism, which is called TAASS, uses a time aver- 
age estimate of the correlation of (\y(i)\ 2 — 1) and (\y(i — 1)| 2 — 1). 
The update rule is 

H(i + 1) = ap,(i) + -yv 2 (i) (8) 

where v(i) = Pv(i - 1) + (1 - P)(\y(i)\ 2 - l)(\y(i - 1)| 2 - 1) and 
< P < 1. The limits on fi(i + 1), ct and 7 are similar to those of the 
MASS algorithm. The exponential weighting parameter /3 governs 
the averaging time constant, namely, the quality of the estimation. 
Previous samples, in stationary environments, contain information 
that is related to determine an accurate measure for the proximity 
of the adaptive beamformer coefficients to the optimal ones. Here, 
p should be close to 1. For non-stationary environments, the time 
averaging window should be small for deleting the deep past data 
and leaving space for the current statistics adaption, so, P < 1. 

There are two objectives for using v(i) here. First, it rejects the 
effect of the uncorrelated noise sequence on the step-size update 1 18 1. 
In the beginning, because of scarcity of transmitters' information, 
the error correlation estimate v 2 (i) is large and so n{i) is large to 
increase the convergence rate and to track the change of input data. As 
it approaches the optimum, v 2 (i) is very small , resulting in a small 
step size for ensuring low misadjustment near optimum. Second, the 
error correlation is generally a good measure of the proximity to the 
optimum. 

C. Computational Complexity and Convergence Analysis 

In this part, both the computational complexity and the conver- 
gence behavior of the proposed mechanisms based on the CCM 
criterion are investigated. 

1) Computational Complexity: The computational complexities of 
the proposed MASS and TAASS mechanisms are analyzed. It is 
well known that the computational complexity of the ASS algorithm 
is a linear monotonic increasing function of the number of sensor 
elements (in AWGN model), i.e., the complexity will increase fol- 
lowing the increase of the number of sensor elements. Therefore, the 
computational complexity becomes very large if the array size is big. 

An important feature of the proposed algorithms is that they only 
require a few fixed number of operations for updating the step size 
compared with that of the ASS method, which is proportional to the 
number of sensor elements. The additional computational complexity 
of the proposed adaptive step size mechanisms is listed in Table HI 

2) Convergence Analysis: Considering the space limitation, we 
just give the range of the step size for convergence. For further 
analysis, we assume that p,(i) varies slowly around its mean value. 
This assumption is approximately true if 7 is small and a closes to 
one, which will be shown in the simulations. Under this condition, 
according to 11161 . we have 

E[f,({)(\y(i)\ 2 - l)y*(i)[x(i) - a H (8 )x(i)a(8 )]] 
= E[fi(i)]E[(\y(i)\ 2 - - a H (e )x(i)a(e )]] 

and 

EMi)(\y(i)\ 2 - l)x(i)x H (i)]w{i) = E[fi(i)]R C CMw(i) (10) 



where Rccm = E[{\y{i)\ 2 - l)x(i)x H (i)] G C mxm . 

Now, the weight vector update equation $5$ of the blind adaptive 
CCM beamformer can be written as 

w{i + 1) = (I- Li(i)(\y(i)\ 2 - l)v(i)x H (i))w{i) (11) 

where v(i) = (I - a(6 )a H (6 ))x(i) G C mxl . 

Define the weight error vector e w (i) and substitute i ll It into the 
expression, we get 

e w (i + 1) = w(i + 1) - id OJ >t 

= (I - v(i)(\y(i)\ 2 - l)v(i)x H (i))e w (i) (12) 

- I 2 _ l)v(i)x H (i)w opt 

Employing {9} and d 10b and taking expectations on both sides of 
J 1 2b . we get 

E[e w (i + 1)] = (I - E[ii{i)]R vx (i))E[e w (i)\ (13) 

where R vx {i) = E[(\y(i)\ 2 ~ l)v(i)x H (i)] = (I — 
a(6 )a H (9 ))RccM and R vx w opt = 0. Therefore, E[w(i)] 
w op t or equivalently, lirrii-joo E[e w (i) — 0] represents the stable 
condition if and only if Yl^oil ~ E[fi(i)]R vx ) —¥ 0. Following the 
idea of the eigenstructure [12] with respect to the correlation matrix 
Rvx, the sufficient condition for d 1 3b to hold implies that 

< EMoo)} < -J- (14) 

"max 

where A^f aa . is the maximum eigenvalue of R vx - 

V. Simulations 

The performances of the proposed MASS and TAASS algorithms 
are compared with other existing algorithms, namely FSS and ASS, 
in terms of output signal-to-interference-plus-noise ratio (SINR). An 
ULA containing m — 16 sensor elements with half-wavelength 
spacing is considered. The noise is spatially and temporally white 
Gaussian noise with power a 2 , = 0.01. For each scenario, K — 1000 
iterations are used to get each simulated curve. In all simulations, the 
desired signal power is a 2 = 1. The BPSK modulation scheme is 
employed to modulate the signals. 

Fig. [2] includes two experiments. Fig. |2ja) shows the output SINR 
of each method versus the number of snapshots, whose total is 1000 
samples. There are five interferers in the system, one interferer with 
4 dB above the desired user's power level, one with the same power 
level of the desired one and three with power 0.5 dB lower than that of 
the desired user. In this environment, the actual spatial signature of the 
signal is known exactly. We set the first element of the initial weight 
vector w(0) equals to the corresponding element of steering vector 
of SOI a(9o). Other parameters are optimized with a — 0.98, 7 = 
10~ 3 , no = 10~ 5 , Umax = 10~ 4 and fi min = 10~ 6 for MASS and 
a = 0.98, P = 0.99,7 = 10" 3 , Mo = 10~ 4 , /i m « = 3 x 10~ 4 and 
fJ-min = 10~ 6 for TAASS. We claim that the parameters for the FSS 
and ASS are tuned in order to minimize the performance, allowing for 
a fair comparison with the proposed algorithms. The results show that 
the proposed algorithms converge faster and have better performances 
than the existing algorithms. The steering vector mismatch scenario is 
shown in Fig.[2jb). We assume that this steering vector error problem 
is caused by look direction mismatch. The assumed DOA of the SOI 
is a constant value 2° away from the actual direction. Compared with 
Fig. [Ha), Fig. |2jb) indicates that the mismatch problem leads to a 
worse performance for all the solutions. The convergence rate of all 
the methods reduces whereas the devised algorithms are more robust 
to this mismatch, especially for the TAASS approach, which reaches 
the steady-state very quickly. 

In Fig. [3] The system starts with 4 interferers, two of which have 
the same power as that of the desired signal and the rest of them 



with the power 0.5 dB lower than the desired one. Two more users 
with one of them 2 dB above the desired user's power level and the 
other 0.5 dB lower than that of the desired user, enter the system at 
1000 symbols. In this condition, the parameters are set to the same 
values as those of the previous experiment except fimax = 3 x 10 -3 
for MASS and \i ma x = 5 x 10 -3 for TAASS due to optimization. 
As can be seen from the figure, SINRs of all the algorithms reduce 
at the same time. It is clear that the performance degradation of 
the proposed ones is much less significant than those of the other 
methods. In addition, MASS and TAASS methods can quickly track 
the change and recover to a steady-state. This figure illustrates that 
the proposed algorithms have faster convergence than the reported 
methods even though they are less complex. The experiment shows 
that the proposed techniques exhibit better performance after an 
abrupt change, in a non-stationary environment where the number 
of users/interferers suddenly changes in the system. 



(a)SNR=20 dB, q=6 users, m=16 (b)SNR=20 dB,q=6 users, m=16,mismatch=2 
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Fig. 2. Output SINR versus number of snapshots for (a) ideal steering vector 
condition (b) steering vector with mismatch. 
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Fig. 3. Output SINR in a scenario where additional interferers suddenly 
enter and/or leave the system. 



VI. Conclusions 
In this paper, two novel adaptive step size mechanisms employing 
SG algorithms have been presented to enhance the performance, 
improve the convergence property and reduce the computational 
load of the previously proposed adaptive methods for blind adap- 
tive beamforming. We considered different scenarios to compare 
the proposed MASS and TAASS algorithms with several existing 
algorithms. Simulation experiments were conducted to investigate the 
output SINR. The performances of our new methods are shown to 
be superior to those of others, both in terms of convergence rate and 
performance under sudden change in the signal environment even 
though they are less complex. A complete convergence analysis of 
the proposed algorithms is under preparation. 
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